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Nonacoustic (low-frequency) instability of combustion of a solid propellant in a jet 
motor is considered. A model of unsteady combustion in motors with channel charges 
is proposed. The model takes into account the change in the temperature distribu- 
tion in combustion products induced by the change in the gas pressure (Mache effect). 
The unsteady gas-release rate and the temperature of combustion products are deter- 
mined on the basis of the Zel’dovich phenomenological approach, taking into account 
the change in the propellant-surface temperature (Novozhilov model) and the flame 
temperature (Gostintsev Sukhanov model). The dependence of the stability limit of 
the solid rocket motor (SRM) on the charge- channel length and pre-nozzle volume of 
the motor chamber is determined. It is shown that the Mache effect leads to a signifi- 
cant constriction (by a factor of 1.5 to 2) of the range of stable combustion parameters 
in the SRM. It is found that an SRM with a channel charge has a narrower range of 
stability in terms of the Zel’dovich parameter k than an SRM with a butt-end charge 
and an identical combustion-chamber volume. For channel charges, the stability limit 
depends mainly on the volume of the pre-nozzle part of the combustion chamber. 



INTRODUCTION 

Low-frequency (nonacoustic) instability of com- 
bustion in SRMs is usually manifested as subsidence 
of propellant burning with its subsequent self-ignition 
(“sneezing” of the motor) or in the form of periodic 
pressure fluctuations with a frequency lower than 10 Hz. 
The knowledge of the stability limit and its dependence 
on the propellant-combustion and combustion-chamber 
characteristics is necessary for motor design. Therefore, 
the mechanism of emergence of instability in SRMs, be- 
ginning from the pioneering works of Zel’dovich and 
Leipunskii [1-3], is of traditional interest for the solid 
propellant combustion theory (see, for example, [4-17]). 

Usually, low-frequency instability is called the L*- 
instability, which is related to the existence of the criti- 
cal value of the reduced length of the combustion cham- 
ber 

L* = Wja, 

where W is the free volume of the combustor and a is 
the nozzle-throat area. 
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L*-instability arises if the free volume of the com- 
bustor is not sufficiently large. Such a situation is typ- 
ical of the beginning of combustion, when combustion 
may be affected by many factors (ignition, erosion, un- 
steady burning, incompleteness of burning, etc.); that 
is why the experimental and theoretical studies are so 
complicated. 

Despite the important applied significance of the 
problem, there is no commonly accepted quantitative 
theory of instability in the literature. The first qual- 
itative explanations were proposed by Pobedonostsev 
(see [12]). He related instability to poor ignition and 
also to flameout in long channels. The first physically 
grounded theory was proposed by Zel’dovich in [1, 3], 
where he showed that oscillations are related to the 
thermal inertia of the combustion wave propagating in 
a solid propellant and to the flow-rate characteristics 
of the combustion chamber. Low-frequency instability 
arises when the characteristic time of thermal relaxation 
of the combustion wave tt = ycjv} {x is the tempera- 
ture diffusivity and u is the burning rate of the propel- 
lant) becomes comparable with the characteristic time 
of pressure relaxation in the combustor 
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Tr = WlARTka, 

where A is the discharge coefficient, R is the gas con- 
stant, and Tb is the temperature of the propellant- 
burning products. 

Most theories of low-frequency instability use the 
assumption of a uniform distribution of the temperature 
of combustion products over the combustor volume (and 
even a constant gas temperature). This assumption 
simplifies the mathematical formulation of the problem 
but ignores the gas-temperature dependence in com- 
pression or expansion on the initial pressure at which 
this gas was formed (Mache effect [18]). It was noted in 
some papers that this fact is principally important (see, 
e.g., [6, 7, 12]). Then et al. [7] were the first to propose 
the theory of L*-instability of SRM, taking into account 
the temperature (entropy) waves in the combustor for 
motors with butt-end charges. An essential assumption 
of the theory of [7] limiting its application is the as- 
sumption of the constant residence time of gas particles 
in the combustor, which is considered as a parameter of 
the problem. The analysis of L*-instability performed 
in the present work takes into account the significant 
variation of the particle residence time in the combus- 
tor in the course of burning and also its dependence on 
the place of gas formation (burning-surface geometry). 
This can be done using the Lagrangian description for 
the gas flow in the combustion chamber [19]. 

The present work is the further development of 
the theory of L*-instability, which was suggested by 
the authors previously for motors with butt-end burn- 
ing charges [17]. The general phenomenological ap- 
proach [1, 4, 5] is used to describe unsteady burning 
of the propellant, which allows one to avoid the known 
limitations [9, 12] arising in choosing a particular model 
of the propellant flame [7]. 



1. MODEL OF THE PROCESS 

FOR BUTT-END BURNING CHARGES 

For convenience of consideration, we first consider a 
motor with a butt-end charge, namely, a charge whose 
burning surface is equidistant from the motor nozzle 
(see Fig. la). We briefly describe the main postulates 
and relations of the theory suggested (see [17] for more 
detail) . 

Our consideration is based on the following as- 
sumptions. 

1 . The gas motion from the burning surface toward 
the nozzle is laminar and one-dimensional; gas portions 
formed at different times do not mix during the entire 
time of their residence in the combustor. 



a b 




Fig. 1. SRM schemes: (a) butt-end charge with 

an equidistant surface from the nozzle; (b) channel 
charge. 

2. The gas pressure in the combustor (p) is identi- 
cal over the volume and changes only in time: p = p(t). 
This assumption makes the principal difference from the 
problem of acoustic instability of combustion, where the 
pressure in the combustion chamber is nonuniform, but 
the mean pressure is constant in time. The change in 
the mean pressure is a second-order effect in the acous- 
tic wave amplitude. In the limiting case of quasi-steady 
disturbances, both formulations of the problem coin- 
cide. 

3. The mass velocity of gas outflow through the 
nozzle (Apa) changes in accordance with the current 
values of pressure p{t) and temperature T{t) of the gas 
particle that has reached the nozzle throat by the time t; 
the nozzle-throat area a is constant in time. 

4. The propellant is homogeneous; the density, heat 
capacity, and thermal conductivity of the condensed 
phase are constant. 

5. The composition and heat capacity of combus- 
tion products are constant. This assumption corre- 
sponds to two limiting cases for the relation between 
the chemical reaction time and the particle residence 
time in the combustion chamber. In the first case, the 
chemical reaction zone width and, hence, the reaction 
time are negligibly small as compared to other spatial 
and time scales of the problem. The composition of 
combustion products and the flame temperature corre- 
spond to conditions of thermodynamic equilibrium. In 
the second case (flameless burning), the relationship be- 
tween the scales is the opposite; therefore, the chemical 
reaction has not enough time to produce a noticeable 
effect on the gas state in the combustion chamber. The 
conditions of realization of these regimes were analyzed 
in [14, 15]. 

Despite the uniform pressure, the temperatures of 
gas portions formed at different times may be signifi- 
cantly different due to the Mache effect. As a result, 
the gas flow along the chamber has nonuniform distri- 
butions of temperature and density. The characteristic 
spatial scale I of temperature nonuniformity in the com- 
bustion chamber, which coincides with the scale of dif- 
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ferent gas portions, is determined by the characteristic 
time of pressure relaxation (r^) and gas velocity in the 
combustor (V): I = VTr- 

Heat exchange between gas particles is ignored due 
to the rather low frequency of pressure oscillations and 
the small residence time of the gas in the combustor [17]. 



Basic Relations 



We use the Lagrangian approach [19] to describe 
the gas motion in the combustor and nozzle. This 
method (in contrast to the traditional Eulerian ap- 
proach) allows us to make the formulation of the prob- 
lem independent of the particular charge geometry. 

Let a gas portion of mass dm at a pressure p{t) be 
formed at a time t near the burning surface during a 
small period dt. The temperature of the gas portion 
formed is equal to the burning temperature Tb(t) deter- 
mined by burning conditions at a given time. Within 
the framework of the phenomenological theory of un- 
steady burning of a solid propellant [1, 4, 5] (assum- 
ing that the gas-phase reactions are completed and the 
gas-phase combustion wave is inertialess), the temper- 
ature Tb is a function of instantaneous values of pres- 
sure p{t) and temperature gradient in the propellant 
near the burning surface 4>(t): Tb = ?b(p, 4)). 

We denote the time of motion of the gas portion to 
the nozzle as r„. Then, at the moment of outflow t + Tn, 
the pressure in the combustion chamber is p{t+Tn), and 
the temperature of this gas portion T{t + Tn) depends 
on Tb(t), p{t), and p{t + Tn)- The residence time r„ in 
an unsteady process is a function of the time t. 

For an adiabatic (in the general case, polytropic) 
process of Lagrangian particle motion, the temperature 
of the portion formed at the time is 

T(t)=Tb(t,)[p(<)Mt,)](^-i)/^, 

where N is the ratio of specific heats. 

Taking into account that the discharge coefficient A 
depends on the gas temperature at the nozzle entrance 
as we obtain 

A{t) = Ah{t - Tn)[p{t - r„)/p(t)](^"^)/2^, (1) 

where Hb(t) is the discharge coefficient at the tempera- 
ture 

The free volume of the combustion chamber occu- 
pied by the gas is 







(2) 



where 7 is the propellant density, p is the gas density, 
S is the burning surface area, and u(t) is the instan- 
taneous burning rate of the propellant. For a given 



combustor volume, formula (2) should be considered as 
an integral equation with respect to p{t) . 

Taking into account that, by definition, the entire 
gas mass formed during the time between t = 0 and 
t — Tn is exhausted through the nozzle by the time t, we 
may write 

t — Tn t 

j dtf. = J A{tA}p{t*)<J dtf,. 



Differentiating this equation with respect to t, we obtain 
^ ^ A{t)p{t)a 
dt u{t — Tn)^S 

System (l)-(3) describes the dynamics of entropy per- 
turbations in the SRM combustion chamber (Mache ef- 
fect) . 



Instability of Steady Burning 

We consider unsteady burning in an SRM with a 
small perturbation of the steady regime: 

uojS = Aopoa. 

It follows from Eq. (2) that the steady value of the 
residence time of the gas in the combustion chamber 
equals the time of pressure relaxation in the SRM com- 
bustor (Tr): 

_o Wpo W __ 

jSuo AoRn^o^ - 

We pass to dimensionless variables 

■Un M , 4> , 

T = —t, V= 1, (f = —~l, 

K Uo 4*0 

P w , Tn o ^s,0 

V= 1, 7^, 

Po t); Ts^o - To 

where Tg is the burning surface temperature and Tq is 
the initial temperature of the propellant. We rewrite 
system (l)-(3) in these variables: 



X ^ v{t) - 

T 

= f {a - ^ ^ ^ ^?7 (th.) - k£y(TH.) -ku(rH.) dr*, (4) 



dip rN+1 . . i / iV - i\ , , 

= - - 2 1 “ - 

- 2 ^t{t - x) - v{t - x) ■ (5) 
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Here x = a is the sensitivity of the flame tem- 

perature to the gas pressure, and v is the relative fluc- 
tuation of the burning rate. 

Substituting perturbations in the form 



? 7 (r) = r]i exp(Hr), v = v\ exp(Hr), 
ip = (fii expiflr), Ip = -ip I exp (fir) 



into (4) and (5) and eliminating 'tpi, we obtain the rela- 
tion between perturbation amplitudes: 



■ 1 
-N 



^X + 



N+1 ( N-1 

2N “ W~ 



)(l-^ exp(-tix)) 

exp(-Hx)^(/?i = 0. 



Vi 

(6) 



Relation (6) closes the system of equations of the theory 
of unsteady burning of the propellant [4, 5] . 

The characteristic equation derived from the condi- 
tion of solvability of the system of equations of unsteady 
burning and from Eq. (6) has the form 



VI + 4ff - 1 r 



2n 






ky + e{6-i^)[l- - exp(-Hx) 



— {k + r — l)y + S — v — 



1 



vmn-h 1 



^exp(-ffx))! =0, (7) 



nx N+1 
^ 2N 

Here 






y = 



k — {Tsft — To) 

1 ^Ts,o 
Ts,o - To dlnp’ 



din uq 
dTn ’ 



r = 



din uo 
dlnp ’ 



dTs,o 
dTo ’ 

5 = vr — yk 



are the characteristics of sensitivity of the steady burn- 
ing rate of the propellant uq and the surface tempera- 
ture Tsp to changes in the initial temperature and pres- 
sure [4, 5]. 



Some results of the numerical solution of Eq. (7) are 
plotted in Fig. 2 (see also [17]). Curve 1 corresponds to 
isothermal combustion products (a = £ = 0) and an 
isothermal process in the SRM combustor (N =1) — 
there is no Mache effect. In the case of a nonisothermal 
polytropic process in the combustor (e.g., for N = 1.2), 
even for isothermal combustion products (a = e = 0), 
the region of stable combustion in the SRM becomes 
narrower (curves 2). Nonisothermal combustion prod- 
ucts additionally destabilize SRM combustion. The sta- 
bility limit for this case is shown by curves 3 (a = 0.25 
and e = 0) and 4 (a = 0.25 and e = 0.25). It should be 
noted that all curves converge rather rapidly (y V 2) to 
a common vertical asymptote (k = k*) corresponding 
to the boundary of stable combustion of the propellant 
at constant pressure [9]: r = (k* — l)^/(k* + 1). 

It follows from the results obtained that the stabil- 
ity limit does not depend on the combustor length and 
is determined only by its volume (more exactly, by the 
reduced length of the combustion chamber L* = W/a). 
At first glance, this result contradicts the model pro- 
posed, since L*-instability may be interpreted (by anal- 
ogy with acoustic instability) as the resonance between 
entropy waves in the combustion chamber and the os- 
cillatory system, which is the burning propellant [9]. 
However, as is shown in [17], this contradiction is ap- 
parent; the reason is that, in butt-end combustion, the 
residence time of the gas in the combustor equals the 
time of pressure relaxation and depends only on the 
combustor volume (and not on its length). 

In the general case (for instance, for channel 
charges), the residence time of the gas in the combus- 
tor depends significantly on the fact at which point of 
the burning surface this or that portion of the gas was 
formed. Therefore, the condition of stability of an SRM 
with a channel charge may be expected to depend on 
the charge-channel length. 



Stability Limits of Steady Combustion 

The boundary of parameters of stable combustion 
of the SRM propellant is determined by solving the 
characteristic equation (7) with imaginary values of the 
frequency H. First, from (7), we can easily obtain the 
critical relation between the propellant parameters at 
which quasi-steady oscillations are formed (H — > 0): 

5 - r(l - a/2) + e(S + p)/2 = 0. 

From here, for e = /u = 0, we obtain the critical relation 
ly = 1 — a/2, which is in agreement with the results 
of [11], where the stability of combustion to quasi-steady 
acoustic oscillations in the motor was studied. 



2. MODEL OF THE PROCESS 
FOR CHANNEL CHARGES 

We consider the scheme of a motor with an ax- 
isymmetric channel charge (see Fig. lb). Propellant 
combustion occurs only over the inner surface of the 
channel: the butt-end faces and the external surface of 
the charge are reinforced. There is a pre-nozzle vol- 
ume between the back end face of the charge and the 
chamber nozzle. 

The following basic assumptions are used. 

1. The pressure in the combustor (including the 
channel) is uniform; the effect of gas-dynamic pressure 
differences on the combustion and exhaustion processes 
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Fig. 2. Stability limit versus the characteristics of steady combustion of a butt-end charge for /i = 0 (a) and 0.2 (b): 
1) isothermal gas in the combustor; 2) isothermal gas release from the burning surface; 3, 4) nonisothermal combustion 
products. 



is ignored (this corresponds to subsonic flows in the 
combustor) . 

2. All portions of the charge surface burn simul- 
taneously, i.e., with an identical unsteady burning rate 
and flame temperature. 

3. Mixing of different portions of combustion prod- 
ucts in the channel, combustor, and nozzle and heat 
transfer between them are negligibly small (the motion 
of each gas portion is adiabatic). 

Each of the above assumptions imposes certain re- 
strictions on the channel length. For example, it fol- 
lows from assumption 1 that the channel length should 
be smaller than the critical value [20] at which sonic 
choking of the channel occurs. Assumption 3 limits the 
channel length by a value corresponding to the begin- 
ning of flow turbulization (the case of turbulent mixing 
of combustion products was considered in [7, 9, 14, 15]). 



Basic Relations 

In the case considered, each elementary portion of 
combustion products is characterized not only by the 
time of its formation t* but also by the point of the 
burning surface, where it was formed. To character- 
ize the position of the point on the burning surface, 
we introduce the variable S, which is numerically equal 
to the surface area enclosed between the channel cross 
section considered and a certain “initial” cross section. 
The initial cross section is chosen in such a way that 



the time r„ is a monotonically increasing function of 
the parameter S: 0 ^ S ^ So, where Sq is the total 
burning surface of the charge, which coincides with the 
channel surface. 

Taking into account that r„ = r„(f. S'), the dis- 
charge coefficient for each gas portion described by for- 
mula (1) is also a function of the parameters t and S: 
A = A(t,S). 

Owing to the assumption of simultaneous combus- 
tion of all points of the charge surface, the quantities u, 
Tb, Ab, and 4> are functions of the time t* only. 

Passing through the nozzle, each portion of com- 
bustion products forms a stream tube with the mini- 
mum cross section da. 

The law of conservation of mass for combustion 
products formed at the section dS of the charge sur- 
face has the form 

t — T„{t,S) t 

dS J dt„ = da J S)p{t^) dt,^. 

0 0 

Differentiating this relation with respect to t, we obtain 
the equation 

r dTr,{t,S)iu{t-Tr,)j^da 

dt \ A{t,S)p{t) ~ dS’ 

Obviously, the sum of instantaneous areas of the 
minimum cross sections of all stream tubes correspond- 
ing to different portions of combustion products is equal 
to the area of the minimum cross section of the nozzle cr. 
Then, integrating Eq. (8), we obtain 
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dTn{t,S) i u{t - r „)7 

dt \ A{t, S)p{t) 



dS = . 



(9) 



Equation (9) is the extension of Eq. (3) to the case of 
channel charges. 

To close the system of equations, it is necessary to 
determine the dependence r„(t, S'). For this purpose, 
we assume that the velocity distribution of combustion 
products over the channel cross section is uniform. This 
assumption is widely used in flow calculations in SRM 
combustors. 

The combustor volume (see Fig. lb) may be di- 
vided into two parts: channel volume and pre-nozzle 
volume Wq. Between an arbitrary cross section of the 
channel S and the “initial” cross section, there is a 
volume W (S) that contains only combustion products 
formed at the section [0, 5] of the burning surface. For 
a given charge geometry, the function W (S') is always 
known. 

By analogy with (2), for the volume VF(S), we ob- 
tain 



s t-Tn(t,S) 



W{S) = 



p(t*)l 7w(t*) 



0 t — Tn{t,S*) 



P{t) J P[p{t^)] 



dt^ dSf (10) 



From the same reasoning, for the total volume of the 
combustion chamber = IF(So) -I- ITo, we can write 



M^(So) + IFo 



So t 




ju{u) 

- P{t) i p[p{Q] 



dt^^dS^,. (11) 



Subtracting Eq. (10) for S = So from Eq. (11), after 
simple transformations, we obtain 

t 



lFo = 



t-T„(t,So) 



p{t*) ]y^ ju{U)So 

P{t)\ p[p{U)] ^ ’ 



which exactly coincides with Eq. (2). 

System (1), (9), (10), (12) describes the dynamics 
of entropy perturbations in the SRM combustor with a 
channel charge. 

Differentiating (10) with respect to S, we obtain 



dW{S) 

dS 



= -S 



- p{t - T„) 1 1/^ 

- P{t) - 



Ju{t-Tn) dTn{t,S) 

p[p{t - T„)] ds ■ ^ 



We confine ourselves to the case dW{S)/dS = const. 
This corresponds to a constant-area channel (without 
the infiation effect [2, 12]). 

It follows from Eqs. (12) and (13) that the following 
formula is valid for the steady regime of combustion: 



Tn(S) = Tn(^o) ~ a — In 

juo So 

VEo 



(14) 



r°(So) = 



Wopo 



= Tr 



7S0M0 AoRTb^oa 

We transform Eq. (13) to a more convenient form 
for further analysis: 

So i-T„(t,So) 





\p{t*)] 


J S dS J 


1 p{t) \ 



■ dt*. 



(15) 



t — Tn{t,S) 



Stability of Steady Combustion 



We consider small perturbations of the parameters 
relative to their steady values: 

i’o = - 1 , ^s(t,S) = [r„(t,S) - r°(S)] — . 

Ty’ ^ 

We define the instrumental constant Xo as Xo = 
TrUg / K. Note that the constant xo depends only on the 
pre-nozzle volume of the SRM combustor. 

We introduce the function 






(16) 



which may be considered as an “instrumental” function 
of the motor, where the parameter 

pouo aSo 

P = a = — Xo 

7>c Wo 

reflects the effect of the channel geometry. 

For small perturbations, Eq. (12) takes the form 
similar to Eq. (8): 



Xo 



- V'o('t) 
Ln 



/[(“- 



N - 1 
N 



p{tA} + eipirAj + v{rAj dr^, (17) 



T-xo 

and Eq. (15) acquires the form 

T-Xs{So) 

iPs{t,So) -iPs{t,S) = J 



a — 



N -1 
N 



-r-XsiS) 

+ ev5*('T*) + r’*(r*) 

From Eq. (9), we can readily obtain 
So 

N+1 



??*(r*) 

dr». 



(18) 



2N 



p4t)So = J v^T - Xs) + 



, ^ 1 ^ dtps{T,S) 

>^{r-Xs) + ^£<P*{t-Xs) ^ 



dS. (19) 
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Substituting perturbations in an exponential form 
and also the perturbation rps(T,S) in the form 
tps{T,S) = ■i/'s,i(<S') exp(rir) into Eq. (18), after simple 
transformations, we obtain 



A,i = xoi^i 



1 ■ 
h . 




N -1 
N 



rji + eipi + vi 



r / S' \ 

l-(^) J exp(-nxo)- (20) 



In deriving this formula, we took into account the ob- 
vious condition ^/>s(r, Sq) = Xo^* or V’s,i(Sq) = XoV’i- 
Substituting Eq. (20) into (17) and (19) and elim- 
inating ipi, we obtain 



■1 

-n 



^Xo + 



N+l 

2N 




IV- 1\ 
N ) 



X 1- 



1 



2(1 -|- /3fl) 
1 



1 - 



2{i + pn) 



exp(-flxo)jj»7i - vi 

exp(-flxo)j£7’i = 0. (21) 



The characteristic equation obtained from the con- 
dition of solvability of equations of unsteady combustion 
and Eq. (21) has the form 



VI + - 1 [ 

2X1 



X 1- 



ky + e{S — v) 
1 



2(1 -|- j3Xl') 



exp(-flxo) -{k + r-l)y 



6 — V — 



\/i + 4f2 + 1 r 



5 — ry + s{5 + n) 



X 1- 



2{i + pn) 



exp(-flxo) 



= 0 , 



(22) 



where 

1 ^ N+l 



N ) 



X 1- 



1 

2(1 + 



exp(-flxo))- 



For /3 = 0, Eq. (22) is transformed into Eq. (7) for an 
SRM with a butt-end charge. It can be easily seen that 
the parameter j3, which reflects the channel effect, enters 
Eq. (22) in the form of the product (3X1. Therefore, in 
the limiting case of quasi-steady perturbations (fl — > 0), 
the critical conditions for channel and butt-end charges 
coincide. In the general case {XI > 0), the critical con- 
dition and the corresponding frequency of oscillations 
depend on the parameter (3. 




Fig. 3. Limits stability for different values of the 
pre-nozzle volume of the combustion chamber with 
a channel charge. 

Stability Limits 
of Steady Combustion 

We consider the effect of the parameters (3 and xo 
on the stability limits in the plane (X)^)- Note that 
the instrumental constant for the free volume of the 
combustion chamber is x = /3 + Xo, where xo and (3 cor- 
respond to the pre-nozzle volume and channel volume 
(see Fig. lb). The results of the numerical solution of 
Eq. (22) at the limit of stability for r = 1/3, v = 2/3, 
= 0, a = 0,25, e = 0, n = 1,2, and different values 
of the parameter xo are shown in Figs. 3-6. The re- 
gion of instability in Fig. 3 is located to the right of the 
curves. It is seen from the figure that the stability limit 
depends not only on the total volume of the combustion 
chamber x but also on the fraction of this volume xo 
that refers to the pre-nozzle volume. 

The curves presented are limited from below by the 
line (/3 = 0, x = Xo)> which corresponds to the limit of 
stability for a butt-end charge. It is seen from Figs. 3 
and 2a that the charge shape has a significant effect on 
low-frequency instability of combustion in SRM. Thus, 
for an identical volume of the combustion chamber, an 
SRM with a channel charge has a significantly narrower 
region of stability in terms of the parameter k than an 
SRM with a butt-end charge, the region of stability be- 
coming narrower as the pre-nozzle volume decreases. 

As it follows from Fig. 3, the stability limit may be 
represented as the dependence k* = k*{x,Xo)- K can 
be easily seen that the dependence on x is weak, and we 
may approximately write k* = k*(xo)- With a sufficient 
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Fig. 4. Dimensionless frequency of oscillations at the 
limit of stability versus the pre-nozzle volume of the 
combustion chamber with a channel charge. 



degree of accuracy, we may assume that the process in 
the SRM is unstable if A: > k*(xo) or xo < X*{k); the 
functions k = k*{xo) and Xo = X*{k) are reciprocal and 
describe the boundary of low-frequency instability for 
butt-end charges. The dependence y = x*{k) is shown 
by the lower curve in Fig. 3. 

According to the data of Fig. 3, we may argue that 
the limit of low-frequency instability for channel charges 
does not depend on the total volume of the combustion 
chamber but is mainly determined by the pre-nozzle 
volume. Therefore, if the combustor volume of an SRM 
with a butt-end charge coincides with the pre-nozzle 
volume of an SRM with a channel charge, their limits of 
low-frequency instability coincide. It follows from here 
that it is necessary to increase the pre-nozzle volume to 
stabilize the process in the SRM with a channel charge, 
since an increase in the channel volume has almost no 
effect on the stability of the process. 

We also note a practically important feature of 
combustion of channel charges: with increasing free vol- 
ume of the combustion chamber due to increasing chan- 
nel volume, the limit of stability in the SRM does not 
approach the limit of stability of propellant combustion 
at constant pressure (which is observed in butt-end com- 
bustion) . 

Figure 4 shows the eigenfrequencies (in a dimen- 
sionless form) for an SRM with a channel charge at the 
limit of stability. Depending on the pre-nozzle volume, 
the eigenfrequency may change severalfold, and the de- 
pendence of the dimensionless frequency on xo is not 
monotonic. 




Fig. 5. Boundaries of the “minimum” stability of an 
SRM with a channel charge for different values of the 
parameter u in the burning rate law. 



As it follows from the results in Fig. 3, the smaller 
the pre-nozzle volume, the narrower the region of sta- 
bility of the SRM. For xo = Oj other conditions being 
equal, the region of stability is minimum. 

Figure 5 shows the boundary of the “minimum” 
region of stability of channel-charge combustion, which 
corresponds to a zero pre-nozzle volume. In calcula- 
tions, we used the same values of parameters as for 
Fig. 3. It is seen from Fig. 5 that an increase in the 
value of ly in the burning rate law makes the region of 
SRM stability narrower in terms of the parameter k. 

Figure 6 illustrates the effect of the parameter iz on 
the frequency of oscillations in an SRM with a channel 
charge at the limit of minimum stability. It is seen 
from Figs. 5 and 6 that the dimensionless frequency 
of oscillations increases with decreasing volume of the 
combustion chamber (curves 1, 2, 4, and 6 in Fig. 4) 
and may be significantly greater than unity if there is 
no pre-nozzle volume (see Fig. 6). 

It is of interest to compare the effect of the gas in 
the combustion chamber with the influence of thermal 
inertia of the gas phase of the combustion wave at con- 
stant pressure [20]. In contrast to the gas phase of the 
combustion wave, the gas in the channel destabilizes the 
combustion process. Nevertheless, the damping action 
of the gas in the pre-nozzle volume of the motor is sim- 
ilar to the effect of thermal inertia of the gas phase of 
the combustion wave. 

We may assume that the Mache effect is one of the 
important reasons for the existence of a critical chan- 
nel length after which channel ignition and subsequent 
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Fig. 6. Dimensionless frequency of oscillations at the 
limit of the “minimum” stability of an SRM with a 
channel charge versus the channel volume and the 
parameter v. 

combustion are accompanied by low-frequency oscilla- 
tions of high amplitude [21]. 

Thus, the model proposed allows one to analyze 
the mechanism of the stabilizing effect of the pre-nozzle 
volume on the combustion process in the SRM, which 
was first observed experimentally by Leipunskii [2]. 

Finally, we note that the proposed model of un- 
steady combustion in the SRM may be useful in analyz- 
ing the Mache effect in motors with more complicated 
shapes of the charge and combustion chamber (for in- 
stance, for multi-cartridge and multi-channel charges, 
and also for channel charges of uniform combustion). 



CONCLUSIONS 

It is shown that the frequency characteristics and 
critical conditions of nonacoustic instability of combus- 
tion in SRMs with channel charges depend on the size 
ratio of the charge channel and pre-nozzle volume of 
the combustion chamber. The nonequilibrium distribu- 
tion of temperature in combustion products along the 
combustor is one of the most important factors that in- 
fluence combustion stability in the SRM, which is com- 
parable with the effect of the dependence of the burning 
rate on the initial temperature. 

For X < 2, the Mache effect leads to an increase 
in the critical value of the instrumental constant y by 
a factor of 1.5 to 2, with an unchanged value of the 
Zel’dovich parameter k. This effect necessitates an in- 



crease in the minimum allowable free volume of the com- 
bustion chamber by a factor of 1.5 to 2, other conditions 
being equal. The main stabilizing effect is exerted by 
the increase in the pre-nozzle volume of the combustion 
chamber. 

With increasing free volume of the combustor due 
to the increase in the channel volume, the combustion- 
stability limit in the motor does not approach the 
combustion-stability limit at constant pressure (which 
is observed in butt-end combustion). 

This work was devoted to the 90th anniversary of 

O. I. Leipunskii. 
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